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Abstract. We calculate the electromagnetic form factor of the pion in quenched lattice QCD. The 
non-perturbatively improved Sheikoleslami-Wohlert lattice action is used together with the &{a) 
improved current. We calculate form factor for pion masses down to m K = 360 MeV. We compare 
the mean square radius for the pion extracted from our form factors to the value obtained from the 
'Bethe Salpeter amplitude' . Using (quenched) chiral perturbation theory, we extrapolate our results 
towards the physical pion mass. 



INTRODUCTION 

The pion, being the lightest and simplest particle in the hadronic spectrum has been stud- 
ied intensively in the past. Using a variety of effective and phenemenological models, 
the properties of the pion have been desribed with varying succes. However, these mod- 
els make assumptions, for example, confinement is put in by hand in contrast to being 
the result of the underlying dynamics. Lattice QCD (LQCD) doesn't have this drawback 
since it is solves QCD directly. 

Using LQCD, global properties of the pion such as the mass and the decay width 
have been calculated to satisfying accuracy. The form factor, which directly reflects the 
internal structure, is clearly an important challenge. The first lattice results were obtained 
by Martinelli and Sachrajda [1]. It was followed by a more detailed study by Draper et 
al. [2], who showed that the form factor could be described by a simple monopole form 
as suggested by vector meson dominance [3]. We extend [4] these studies by adopting 
improved lattice techniques [10-14], which means that we include extra operators in 
order to systematically eliminate all the G[d) discretisation errors. 

Some aspects of the pion structure have been obtained [5-9] using 'the Bethe-Salpeter 
method'. We also use this approach and compare its predictions to the results of our 
direct calculation of the pion form factor. 

Finally we study a chiral extrapolation to reach the physical limit. 



THE METHOD 

To extract the form factor we calculate the two- and three point functions of the pion, 
analogously to [2]. The two point function is given by 

G 2 (?,p) = ^(^(?,x)^( ,0)) e i P- x , (1) 



where <j>* is an operator creating a state with the quantum numbers of the pion. By 
varying the interquark distance at the sink, tf, we can improve the overlap with the 
physical pion and obtain information on the 'Bethe-Salpeter amplitudes' . The three point 
function is calculated as 

G 3 (^;p / ,p i ) = £ (*(*/)*(*) ^(O))*"'^-*)-'** (2) 

X/ ,X 

with 74 the fourth component of the current, inserted at time t. Since the local current 

jfa) = 9(x) 7m (3) 

is not conserved on the lattice, one can construct the Noether current belonging to our 
action 

j£ = k (v/(x)(l - r M )C/ M (x)y(x+ A) - W(x+fi)(l + y M )t/t . (4) 

This current however, still has 6(a) corrections for Q 2 > 0. A conserved and improved 
current can be constructed [12-14] 

jji = z v {ill (x) +ac v d v T^ v } , (5) 

with 

Tn V = qr(x) i (7 MV Y(x) , (6) 
Zy = Z v (1 + aby m q ) . 

The bare-quark mass is defined as m q = ^- (l/?c — 1/jc c ), where K c is the kappa value in 
the chiral limit and a denotes the lattice spacing. For our simulation, K c = 0. 13525 [15]. 
Comparison of the currents will give us information on the importance of improvement. 

DETAILS OF THE CALCULATION 

On a 2A 3 x32 lattice, we generated ^(100) quenched gluon configurations at /3 = 6.0. 
Subsequently, we calculated non-perturbatively improved (csw = 1-796 [10]) quark 
propagators for 5 different values of the hopping parameter, K. These propagators were 
then combined to two- and three point functions for the pion with masses ranging from 
360 to 970 MeV. 1 For the improved current, we use the parameters Zy, by and cy as 
determined by Bhattacharya et al. [17]. To extract the significant parameters from our 
numerical data, we use the following parametrisations. For the two point function, we 
have 

G2(t,p) = t ^(p)Z«(p){ e -^ + e - £ p^-)} , (7) 

n=0 



The lattice spacing a = 0.105 fm is taken from [16]. 



including the contribution of the ground state (n=0) and a first excited one (n=l). The 
Z% denote the matrix elements, 

Z£(p) = |(Q|^Kp)| 2 , (8) 

and are related to the 'Bethe-Salpeter amplitude', 4>(7?) = y/Z%(0) /Zj(0). 

E® and £p are the energies of ground and excited state, respectively; R denotes the 
interquark distance. The three-point function is parametrised as 

G3(f/,f;p/,p0 =F(Q 2 ) v /z«(p7) Z^y^f 
+ | v /zi(p / )Z (p i )<l, P/ |j M (0)|0,p / >r^/ ( ^ ) - £p V +(1-0)} . (9) 

The production of pion pairs, 'wrap around effects' being due to the propagation of 
states beyond tf and an elastic contribution from the excited state were ignored in our 
parametrisation after being estimated negligable. 

All parameters in the two- and three-point functions - energies E, Z-factors and the 
form factor F(Q 2 ) - were fit simultaneously to the data per configuration. For the three- 
point function, we let the current insertion time t vary from to tf. The value for 
the parameters and their error in these simultaneous fits are obtained through a single 
elimination jackknife procedure. 



RESULTS 

As a byproduct of our simulations we also obtain pion masses for the 5 different K- 
values. They agree with the literature. We also checked the energie-momentum relation 
and up to the energies involved we found that a continuum relation provides the best 
description. These non-trivial tests indicate that our simulations are done correctly. 

Using the different currents, we extract the form factor for our five K-values. Com- 
parison of the results for the Noether current and the improved current yields that the 
effect of improvement can be as large as 25% for the highest momenta considered. The 
improved results are shown in Fig.l. The high accuracy for the datapoint at Q 2 = is 
due to the fact that we satisfy the Ward-Takahashi identity to 1 ppm. 

As in the previous study [2] of the pion form factor we compare our results to a 
monopole form factor 

F(Q 2 ) = {i + ^}-\ (10) 

which is suggested by the vector meson dominance ansatz. Fitting our data to this model, 
we extract a vector meson mass which is within 5% of the corresponding rho mass on 
the lattice [15]. From the behaviour of the form factor at low Q 2 , we can extract the 
mean-square charge radius of the pion, 
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FIGURE 1. Form factors as a function of Q 2 for the five pion masses. Curves: monopolefits to the data. 



where in the last step we assume Eq. 10 and use the fitted parameter My. The results 
are shown in Fig. 2 as a function of the pion mass. Previously, the 'Bethe-Salpeter- 
amplitude' 4>(7?) has been used to obtain estimates of the charge radius, 

V / BS ' 4 f d 3r<t> 2 (\r\) ' K } 

The results based on this procedure are also shown in Fig. 2. As can be seen these val- 
ues are much lower than the actual values obtained from the form factor. Moreover, the 
Bethe-Salpeter results are almost mass independent, in accordance with the observations 
of Refs. [5-7, 9]. This is a known[7] deficit of the approach, which Fig. 2 makes quan- 
titative. We extrapolate our results obtained with Eq. 1 1 using three different parametri- 
sations. First, we try chiral perturbation theorie (#pt). At one-loop order the prediction 
for the radius [18] is 

/ 9\one-loop 9 /i->n 

v )*ft =ci+c 2 \ogm l 7t (13) 

In quenched ^pt the rms is constant to this order. There are however indications that a 
mass dependence appears at higher order [19]. For our masses we restrict ourselves to 
a term linear in m\ [20]. Lastly, we have also used the VMD prescription and assume a 
linear dependence of My on m\. These three expectations are also plotted in Fig. 2. 
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FIGURE 2. Radius of the pion as obtained from Bethe-Salpeter amplitudes, Eq. 12, and from the form 
factor, Eq. 11. Also shown are three different parametrizations of (r 2 ) (see text). 
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